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Environment-Mediated Quantum State Transfer 
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We propose a scheme for quantum state transfer(QST) between two qubits which is based on 
their individual interaction with a common boson environment. The corresponding single mode 
spin-boson Hamiltonian is solved by mapping it onto a wave propagation problem in a semi-infinite 
ladder and the fidelity is obtained. High fidelity occurs when the qubits are equally coupled to the 
boson while the fidelity becomes smaller for nonsymmetric couplings. The complete phase diagram 
for such an arbitrary QST mediated by bosons is discussed. 
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I. INTRODUCTION 



One of the most interesting problems in the area of 
quantum information is how to transfer a quantum state 
from one location to another. For example, QST from A 
to B can be done via quantum teleportationi if a prior 
connection between the remote places is established, by 
letting A and B have one from two auxiliary entan- 
gled qubits. A more direct method should involve fly- 
ing qubits, such as photons in optical fibers, which can 
send directly quantum information between two distant 
locations of a quantum computer^. However, the latter 
approach would be very difficult to realize experimentally 
since it requires an interface between the optical system 
and the hardware where a quantum computation takes 
place. A very successful QST was pioneered in^ which 
includes the locations A and B and also the quantum 
transfer ground into the same system. In this case the 
quantum information of the qubits at one end of the chain 
propagate via the interaction between the components 
of a permanently coupled physical system or a quantum 
graph. A perfect or nearly perfect QST occurs between 
two local spins in quantum spin chains and networks, 
at least for short distances^^. The protocol for such 
quantum communication relies on the fact that the aux- 
iliary device which plays the role of the quantum channel 
or quantum bus is the medium itself. For example, the 
exchange interactions of a quantum spin chain can allow 
transfer of quantum information between qubits which 
belong to the first and last local spins while the rest of 
the chain acts as the quantum channel. In this case the 
exchange of information is achieved via magnon elemen- 
tary excitations. Physical devices for such QST could 
be built, at least in principle, and since no external con- 
trol is required they can overcome possible decoherence 
mechanisms. 

In this paper we show that a boson environment could 
be used to transfer efficiently a quantum state by acting 
as a quantum channel. It is known 7 - that entanglement 
can be introduced between two qubits if both are inde- 
pendently coupled to a common heat bath with many 



degrees of freedom. We shall show that even the sim- 
plest possible boson environment which consists of one 
mode can also provide an efficient QST mechanism. For 
this purpose a spin-boson Hamiltonian is introduced^ 
known for many applications in physics and chemistry. A 
related spin-boson model allowed Caldeira and LeggettiS 
to study decoherence via dissipation through a weak cou- 
pling of the spin to many bosons, representing a univer- 
sal realization of a physical environment. Due to weak 
spin-boson interaction the excitations within the boson 
heat bath could be ignored and the problem was solved, 
leading to decoherence^. Our spin-boson model can be 
regarded as an extension of 7 - where two qubits coupled to 
a common heat bath become entangled with each other. 
We show that despite the absence of a direct interac- 
tion between them their coupling to a simple boson en- 
vironment mediates an efficient QST. Environment me- 
diated quantum control for a related multi-mode system 
has been performed inii. 



The proposed spin-boson model allows high fidelity 
QST between two distant locations by choosing suitable 
parameters. In order to make the problem tractable we 
chose the simplest possible quantum channel which con- 
sists of a single-mode boson environment. This is the first 
approximation to a full multi-mode Hamiltonian consid- 
ered ir£ by replacing the coupling to many modes by a 
coupling to an effective boson. Our study proceeds as 
follows: In chapter II the proposed spin-boson model is 
introduced with a double two-level system Hamiltonian 
coupled to a single boson. In chapter III a formula is 
derived for the fidelity of a QST which is obtained by 
mapping the system onto a wave-propagation problem 
in a semi-infinite ladder. The results of our calculations 
with the display of the corresponding phase-diagram and 
a discussion about the efficiency of the scheme follow in 
chapter IV. Finally, in chapter V we discuss possible ex- 
tensions and applications. 



2 




E 

FIG. 1: The proposed protocol for QST between two qubits A 
and B represented by two-level systems which interact with a 
common bosonic environment E acting as a quantum transfer 
channel. 

II. MODEL AND AVERAGE FIDELITY 

The studied system is illustrated in Fig. [T] The qubits 
A and B are not directly coupled with each other but are 
connected via an auxiliary boson environment E both 
having nonzero interaction with E. The qubits in A and 
B can be represented by two local spins and E acts as 
the quantum channel. Of course, if E is replaced by a 
quantum spin chain the model reduces to that studied 
The Hamiltonian is given by the sum 

H = u A a A + uj%a% + uitfb 

+ \ A (b + tf)a A + \ B (b + tf)a x B , 

with the qubits in A and B modeled by two-level systems 
of separations u> A , w A , the quantum channel described 
by a single-boson mode environment of frequency u> and 
nonzero linear couplings X A and Xb exist between the 
qubits and the boson channel E, with a x l y l z the corre- 
sponding Pauli matrices. Note the similarity of H to a 
multi-mode model used to study entanglement between 
the qubits in quantum control theory-ii. The main dif- 
ferences between the present study and^iii lies in the 
number of modes and the presence or not of couplings 
between the qubits and the quantum channel. We con- 
sider nonzero spin-boson couplings X A and Xb since they 
are expected to be comparable to the two-level separa- 
tions uj a and tu B . 

The single-mode Hamiltonian H although simple 
enough it cannot be solved exactly. The Hilbert space 
consists of a direct product of three parts with basis 
states \t]a, Vb, m), where, rj A / B = 0,1 label the qubits 
and m — 0,1, 2, 3, ... is the single phonon excitation num- 
ber of the states in the quantum channel. The QST in 
this system can be studied similarly to that in a spin 
network^. Suppose that at time t — an unknown state 
|V>a) = cos(0/2)|O) + e** sin{6 /2)\l) with parameters 9, 
<fi, is generated at qubit A and has to be transferred to B. 
We also initialize the state of the qubit B to |0) and the 
state of the quantum channel E to its lowest boson state 



|0). The initial state of the whole system is \ipA,0, 0) 
which is separable. When evolution takes place the final 
state at time t in general becomes a non-separable mixed 
state. The measurement of the state of qubit B is de- 
scribed by its reduced density matrix and both efficiency 
and quality of the quantum communication is obtained 
by evaluation of the corresponding fidelity 3 . 

The fidelity is usually computed by taking average 
over all pure input states \4>a) in its corresponding Bloch 
sphere 

(Ht)) = ^J dn{i> A \PB(t)\i>A), 

where the state of A to be transferred is \iPa), Ps(i) 
is the reduced density matrix of the qubit B at time t 
and the average is over all initial \iPa)- If we let the 
system evolve for a period of time t m , one can find the 
maximum average fidelity {T) m from the time taken for 
the average fidelity to reach its first peak corresponding 
to the maximum fidelity. The peak time t p is the second 
important quantity which can characterize a quantum 
channel, the first being the average fidelity {T) m . High 
fidelity implies better quantum channel for QST while 
shorter time to reach the peak means faster QST. If (J-) m 
becomes exactly unity we have perfect QST— with the 
quantum state transferred from A to B without any loss 
of quantum information. 

The reduced density matrix for qubit B can be written 

p B = Tr A ,E[pt] 

by tracing out A and E of the evolved total density matrix 
pt = U(t,0)p U(0,t), with initial value po = Pa®Pb®Pe 
and time evolution operator U(t,Q) — e lHt , U = 1. This 
allows to calculate the average fidelity for any time t, 
which we shall simply call it fidelity from now on. As it 
stands this formula is rather complicated to perform an 
analytic evaluation. In the next chapter the problem is 
mapped onto an equivalent wave propagation involving 
two ladders and the corresponding fidelity is written as 
a function of waves propagating in these ladders. 

III. WAVE PROPAGATION 

A parity symmetry present in H simplifies the 
Hamiltonian 9 making it block-diagonal in a suitable two 
qubit Bell states basis 

|*±,m) = -=(|Q0m) ± |llm)) 
v2 

|$±,m) = -=(\Qlm}± \10m)). 



The states split into two having zero matrix elements 
between each other and the block-diagonal Hamiltonian 
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FIG. 2: The equivalent wave propagation in "+" and "-" 
ladders. The nodes denote the basis states and the lines the 
hoppings, the green lines denote ^Jm(\A + As), the blue lines 
^/m{\A — As), the red lines u>% the black lines lj% 

and the on-site energies are com. 



(a) 



matrix is illustrated via two decoupled ladders in Fig. 
[5J The states are represented by nodes and hoppings 
between the nodes by the connecting lines. Note that the 
ladders of Fig. 2 are rather similar to each other, their 
only difference being the ordering of red and black lines. 
This becomes very helpful for our calculation given in the 
Appendix where the computation is shown to simplify in 
the chosen basis. 

The obtained formula for the fidelity can be given in 
the form 

^ = Ta E( Tr K^™] + Tr[BlB m ] + Tr[ClC m ]) 



with 



T z2 



= f m ( + ,t)+* Z fm(~,t) 

= fm(+,t) + (-ia y )f m {-,t)a x 
+ ^r.U-,t) + a z2 ,f m (-,t)<j x 



1 



-1 



where a +, ~' z are the Pauli matrices and f m (±,t) is the 
propagator in the ladders shown. In the notation used, 
e.g. /3(+,2) means the propagator from to = to slice 
to = 3 at time t — 2 in the ladder with parity " +" . This 
gives the fidelity of QST written as a linear combination 
of the propagators in each of the two ladders. 

Since both ladders are semi-infinite the corresponding 
Hilbert space must be truncated at a maximum phonon 
number to. In order to approximate propagation for very 
long times long ladders with large maximum m are re- 
quired. However, a careful study of the formula shows 
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FIG. 3: The convergence of the average fidelity J-(t) as a 
function of the maximum phonon number m: (a) symmetric 
case with Xa = As = 1.0 and lua — ujb — us where the 
parameters in parentheses on the right of the figure are (t, u>s), 
and (b) non-symmetric case where the parameters displayed 
in parentheses are (Xa, As, wa, ujb)- Time t is set to 32000. 



that the fidelity simply arises from the difference between 
propagators in the two ladders. For example, for lua = 
or lob = the two ladders are exactly the same and the 
fidelity becomes precisely zero. Since their structure is 
rather similar, except for the ordering of lines, if a wave 
reaches very far from the origin in one of them a very 
small difference between the two propagators is expected 
with no contribution to fidelity. Therefore, accurate com- 
putations of fidelity do not require very long ladders and 
reasonable maximum to suffices, as seen in Fig. 3. 

The accuracy of the computed results is shown in Fig. 



3(a) by plotting the fidelity as a function of the maximum 
phonon number to for the symmetric case with Xa = 
As = As, to a — <jJb — ws and in Fig. |3(b)| for the 
non-symmetric case. The fidelity is shown to converge 
very rapidly for maximum phonon numbers to = 40 or 
50 which permit to use reasonable coupling strengths. 
The convergence does not depend on time t m and is also 
rather insensitive to loa/b since it mostly depends on the 
couplings Xa/b- F° r example, the numerical results for 
Xs = 1.0 and As = 2.0 required only to = 40 to 50 and 
more that to = 100, respectively. In our computations 
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suitable maximum m was chosen according to the values 
°f ^a/b and the convergence was checked by varying m. 
For couplings A = 0.0 to 2 and u> = 0.0 to 80 a maximum 
phonon number m between 50 to 110 was sufficient. 

IV. RESULTS AND DISCUSSION 

The quality of QST is determined by the maximum of 
the average fidelity {T) m in the time period from to t m 
and the time for occurrence of the first peak t p when the 
system reaches its maximum. Higher fidelity means more 
faithful state transfer while shorter peak time t p implies 
faster state transfer. The parameters wa, ujb, Xa, Xb 
are taken in units of uj — 1 while the maximum fidelity 
and the first peak time are obtained in the time interval 
[0,t m = 32 x 10 3 ]. 




(a)The phase diagram of the maximum fidelity as a function of equal 
couplings Xg and equal qubit separations ajg. Three regions can be 
distinguished as explained in the text. 



A. Phase Diagram for Symmetric Couplings 



The phase diagrams of the maximum fidelity and the 
first peak time are shown in Figs. 4(a), 4(b) as a function 
of the two parameters = u>a = ujb and As = Xa — 
Xb- They can be divided into the following three regions: 

Region I: a weak coupling region which lies in the up- 
per left corner of Fig. 4(a) where ujs » As- In this 
case the corresponding first peak time t p shown in Fig. 
4(b) is large equal to the upper bound of the chosen time 
interval t m . In other words, the fidelity never reaches its 
maximum within the adopted evolution time. This indi- 
cates that probably a higher fidelity might occur for even 
longer times so that we can call this a " slow region" . We 
may conclude that a good state transfer is impossible in 
this region because of the long times t p . 

Region II: lies in the lower part of the figure, which is 
too small to be seen in Fig. 4(a) and this plot is magnified 
in Fig. 4(c). In this region lus and As are of the same 
order of magnitude so that the fidelity is again low but 
for a different reason than that of region I. The first peak 
time in this case from Fig. 4(b) is less than t m and the 
QST is affected by increasing wg. For zero no QST 
is possible while it becomes better when increasing the 
qubit-environment coupling As. 

Region III: The rest of Fig. 4(a). One can see that 
in the majority of this region high fidelity occurs with 
the first peak time mostly being less than 5 x 10 3 . This 
region corresponds to a two-valley Hamiltonian and the 
system behaves as a good quantum channel. 



B. Phase diagram for non-Symmetric Couplings 
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(b)The first peak time of the QST as a function of the equal 
couplings Xg and the qubit separations u>g. This picture also has the 
three regions mentioned in the text. 
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(c)The Region II of Fig. 4(a) shown in more detail. 
FIG. 4: Phase Diagram for the fidelity of QST 



We have also considered the non-symmetric case where 
the two couplings and the two frequencies are not equal. 
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FIG. 5: The maximum of the average fidelity (J-)m = 0.998, 
a value close to a perfect QST, becomes lower for deviations 
from equal frequencies io a = wb = 20.0 and equal couplings 
\a = A_g = 0.8, with Suj — lob — lja and SX = As — A,4, 
respectively. The {3~)m is obtained in the region [0, 33 x 10 3 ]. 



The influence of a deviation from equal couplings is stud- 
ied by choosing As = 0.8, los = 20.0 with the correspond- 
ing point of the symmetric phase diagram belonging to 
region III having very high fidelity equal to 0.998. A 
small deviation Suj in wg with ujb = uja + Sw is shown 
in Fig. 5(a) to influence dramatically the QST, which 
is extremely sensitive even for deviations of the order of 
10~ 4 . The asymmetry in the coupling constants is shown 
in Fig. 5(b) to have a much smaller effect. 



CONCLUSIONS 



Although the role of an environment is usually that of 
causing decoherence for a quantum system the presence 
of entanglement between the system and the environment 
also signals the possibility that quantum information can 
be transferred via the environment. We suggest a QST 
between two qubits via a coupling to a common boson 
medium which acts as the quantum channel. We have 
derived a formula for the corresponding fidelity of the 



state transfer by mapping this problem into a wave prop- 
agation, which is much easier to understand and solve. 
For symmetric couplings and frequency separation case 
high fidelity QST between the two qubits is obtained for a 
wide range of parameters. We show that small deviations 
from this symmetry can dramatically lower the QST. 

Questions for further study are: (i) possible extensions 
of the present scheme to include a multimode boson en- 
vironment since our results can cover only approximately 
the multimode case, (ii) connections of QST to wave 
propagation in media also in the presence of disorder 
which can also give ballistic, chaotic and even localized 
states (in the latter case QST is impossible) and (iii) pos- 
sible realization of an experiment where QST mediated 
by bosons can occur, for example, between two quantum 
dots coupled to the appropriate phonon environment of 
a nanostructure. 
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VII. APPENDIX: DERIVATION OF THE 
FORMULA OF FIDELITY 



The average fidelity 



over \ipA) becomes 
Wa) - 



1 

4tt 



dn(ip A \PB(tMA) 



cc»(-)|Q)+e*«n(^|l) 
dsin{9)d9 / 



l - J dQ.. = ^ dsin{0)d0 



4tt 

The reduced density matrix pB(t) can be calculated via 

Pb = Tr A ,E[pt] 
Pt = U(t,0)p U(0,t) 



U(t,0) 



iHl 



where, the partial trace over the degrees of freedom for 
qubit A and the quantum channel E is taken. H is the 
Hamiltonian for the system A® E and U (t, 0) is the 
corresponding time evolution operator. 

To simplify the formula first we have calculated the 
integral. It is convenient for us to choose coherent vector 
representation^ to express the density matrix. 

PB{t) = 



PA 



(1+pA-ff) 
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an assuming the relation between two coherent vectors 

p B {t) =T{t)-p A + f (t) 
we can carry out the integral 

{T(t)) = \\l + \Tr{T(t))]- 

We need to calculate the matrix T(t), e.g., to express 
the final state of qubit B as a function of initial state of 
qubit A 

p B (t) =Tr AiE [U{t,0)p{0)U(0,t)} 

where p(0) is the initial state of the whole system (A ® 
B ® E), it is separable so that 

p(0) = P a(0) ® Pb (0) ® P e(0) 
Pb(0) = |0)(0| 
Pe(0) = |0)(0|. 



By inserting \riA,riB,m) into these formulae we find 

PB{VB,VB,t) = J B(VB,VB,VA,V' A,t)PA{VA,v' ' A,ty 

VA-.V A 

Jb(vb,Vb,Va,v' A,t) = 



J ^atib , VaVb . VaO, r) A Q, 0) 

VA 



J(.VaVb,VaVb^VaVb,V aV' b,0) = 
^tos, 0|t/(0, Ol^r/s, 0)<^f?B, 0|tf(t, 0)|r/> B , 0), 

where r^/B = 0/1, m = 0, 1, 2, 3, ... 

The matrix element between ps(t) and ^(0) is related 
by the function Jb 

/ T x (01) + T X (10) i[T x (10) -T x (01)} T X (00)-T X (11)\ 
T= T»(01)+T"(10) i[T»(10) - T«(01)j T»(00) - T» (11) 
\ T*(01) + T Z (10) i[T z (10) -T z (01)} T z (00) - T z (ll) J 

where, 

T x (m) = j s (oi,W) + ^s(io,W) 

T»( W ') = i[Js(10,W)-^s(01,W)] 
T z (w') = Js(00,W)-^s(H,W) 

By going into Bell basis the final expression for the fi- 
delity is obtained. 
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